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ABSTRACT 



Gamma-Ray Bursts (GRBs) are traditionally divided to long and short ac- 
cording to their durations 2 sec). It was generally believed that this reflects 
a different physical origin: Collapsars (long) and non-Collapsars (short). We 
have recently shown that the duration distribution of Collapsars is flat, namely 
independent of the duration, at short durations. Using this model for the dis- 
tribution of Collapsars we determine the duration distribution of non-Collapsars 
and estimate the probability that a burst with a given duration (and hardness) 
is a Collapsar or not. We find that this probability depends strongly on the 
spectral window of the observing detector. While the commonly used limit of 2 
sec is conservative and suitable for BATSE bursts, 40% of Swift 's bursts shorter 
than 2 sec are Collapsars and division 0.8 sec is more suitable for Swift . We 
find that the duration overlap of the two populations is very large. On the one 
hand there is a non-negligible fraction of non-Collapsars longer than 10 sec, while 
on the other hand even bursts shorter than 0.5 sec in the Swift sample have a 
non-negligible probability to be Collapsars. Our results enable the construction 
of non-Collapsar samples while controlling the Collapsar contamination. They 
also highlight that no firm conclusions can be drawn based on a single burst and 
they have numerous implications concerning previous studies of non-Collapsar 
properties that were based on the current significantly contaminated Swift sam- 
ples of localized short GRBs. Specifically: (i) all known short bursts with z > 1 
are most likely Collapsars, (ii) the only short burst with a clear jet break is most 
likely a Collapsar, indicating our lack of knowledge concerning non-Collapsar 
beaming (iii) the existence of non-Collapsars with durations up to 10 sec impose 
new challenges to non-Collapsar models. 



Kouveliotou et al 



1. Introduction 

(119931 ) have shown that gamma ray bursts (GRBs) can be divided 
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to two groups according to their observed duration. Long bursts (LGRBs) with observed 
durations T 90 > 2 sec and short ones (SGRBs) with T 90 < 2 sec. They have also found that 
SGRBs are harder on average than LGRBs, supporting further the possibility that the two 
populations arise from different physical sources. Later on afterglow observations enabled 
the localizations of GRBs and identifications of their hosts. These observations supported 
further the different sources hypothesis. Hosts of LGRBs have a large star formation rate 
while SGRB hosts include both star forming and non-star forming galaxies. The position 
distribution of LGRBs within th eir host, towards the center and within high star form- 



ing regions (IFruchter et al.ll2006l ). differs from the position distribution of SGRBs within 



(Barthelmv et al. 


2005; 


Fox et al. 


2005; 


Gehrels et al. 


2005; 


Nakar 


2007; 


Bereer 


2009) 



These observations have led to the realization that GRBs have two different progenitors 
and they are generated by at least two different mechanisms^. The association of LGRBs with 
star forming regions and in sev eral cases with type Ic SNe suggest that they involve stellar 
collapse. The Collapsar model (jPaczynskilll998l ; iMacFadyen fc Woosleylll999l ) suggests that 
a central engine within the collapsing star (powered most likely by an accretion disk onto 
the newly formed compact object or by a magnetar) powers a jet that penetrates the stellar 
envelope and produces the observed gamma-rays once it is outside the star. SGRBs are 
typically weaker and are observed at lower distances. They are more numerous (locally), 
but being harder to detect there are less observed SGRBs than long one. Those SGRBs 
with identified locations are associated with a wide range of stellar populatio n ages. Their 



prope rties are consistent with those expected from bi nary neutron star mergers (jEichler et al. 



19891 ) . although the exact origin is still uncertain (see lNakarll2007l . for a recent review). Since 
the origin of this group is still uncertain we will denote them simply as non-Collapsars. 

It is commonly implicitly assumed that there is one to one correspondence between 
the observed groups of LGRBs and SGRBs and the astrophysical groups of Collapsars and 
non-Collapsars. However, a quick inspection of the duration distribution (fig. [1]) suggests 
that this is not the case and there is a significant overlap between the two groups of long 
and short GRBs: there are SGRBs of Collapsar origin and vice versa. Apart from a slight 
difference in the average hardness, all other high energy emission properties of LGRBs and 
SGRB s are remark ably similar. This makes it difficult to identify the origin of any individual 
burst (lNakarll20071 ) and lacking a better criteria the original division according to T90 2 sec 
is widely used. However this criteria was established for a specific detector (BATSE) with a 
specific observational window. SGRBs are typically harder than long ones and as such they 



1 We have recentl y shown that low lumino sity GRBs are a third group, generated by a different mechanism 
than reg ular GRBs (jBromberg et al.lboilrj ). 
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are more difficult to detect by softer detectors like Swift /BAT than by BATSE. Indeed, 
Swift's short/long detection rate is 1/10 vs. BATSE's 1/3 (when the criterion T go ^ 2 is 
used). This suggests that Swift s division line between the two groups might be at a shorter 
duration, as indeed is seen in a visual inspection of Swifts duration distribution (see fig. [T]). 



Zhang et al.1 (120091 ) suggested to classify individual GRBs using various subsets of prop- 
erties, e.g. spectral lag, peak energy, etc. These subsets of properties are selected phe- 
nomenologically, and are not based on any physical model. The main problem of those 
classification criteria, which are based on the high-energy emission alone, is that there is a 
significant overlap, which cannot be quantified, between Collapsars and non-Collapsars in 
all of them. As a result, the quality of the classification of any of these phenomenological 
methods cannot be quantified and it is therefore impossible to estimate the fraction of mis- 
classified GRBs. This poses a major problem in using such a method especially since the 
sample of GRBs with "good data" (afterglow detection, good localization, redshift measure- 
ments etc.), is small and very sensitive to misclassification. Other attempts used a statistical 
approach and tried to evaluate the overlap between the two populations by fitting the distri- 
bution of GRB s with two underlying distributions. In this case two lognormal distributions 



(IHorvathl 120021 ; iLevesque et al.ll2010l ). The quality of such classification schemes depends 
entirely on similarity of the true distribution of the two populations to the arbitrarily chosen 
fitted distribution. Such approach can be trusted only if we know, for example based on 
physical arg uments, what is the und erling; distribution of at least one of the two populations. 
Recently, in Bromberg et al. (j2011a ) we have shown, based on generic physical properties of 



the Collapsar model, that at short durations the Collapsar distribution is flat. Namely the 
number of Collapsars per unit duration at short durations is independent of the duration. 
Building on this result we estimate here the probability that a GRB with a given duration 
is a Collapsar or not. Not surprisingly this probability depends on the detector and we cal- 
culate it for the three major GRB detectors: BATSE, Swift and Fermi GBM. An improved 
version of this method is obtained by adding a hardness dependence. We present a refined 
probability distribution that is based on both the duration and the hardness. Needless to 
say our method is statistical in nature. We cannot determine whether a specific burst is a 
Collapsar or not, but we can give a probability estimate for this question. 

We begin, in section 2, with a discussion of the Collapsars' duration distribution and 
an analysis of the observed duration distribution of GRBs . In section 3 we calculate the 
probability that an observed GRB is a non-Collapsar. Our results imply that short duration 
Collapsars have been wrongly classified as non-Collapsars, mostly in Swift sample, and this 
has lead to potential misinterpretation of some of the observed data. In section 4 we discuss 
the consistency of our findings with some of the recent studies and the implications on the 
inferred properties of non-Collapsars. We summarize our results and their implications in 
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Fig. 1. — The log(T 90 ) double humped durations distributions, dN / d\og(T 90 ) , of BATSE 
(black) and Swift (red), binned into equally spaced logarithmic bins. Bins with less than 
5 events are merged with their neighbors to reduce statistical errors. The minima in the 
distributions occur at T90 = 2.4 ± 0.4 sec in the BATSE distribution and at T90 = 1.2 ± 0.2 
sec in Swift . 
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section 5. We provide a table of the probabilities for each one of the observed Swift bursts 
with T 90 < 2 sec to be a non-Collapsars in an Appendix. 



The observed GRB duration distribution 



Within the Collapsar model a GR B can only be produced after the jet has emerged from 
the surface of the collapsing star. We (IBromberg et al.ll2011al ) have recently shown that this 
leaves a distinctive mark on the observed duration distribution: it is flat at durations shorter 
than the typical breakout time of the jet from the star (about a few dozen seconds modulo 
the redshift ). In a nutshell, this result arises from a simple fact. The burst duration is 
the difference between two quantities: the engine operating time and the jet breakout time. 
Under quite general conditions the resulting distribution is flat at durations that are shorter 
than the typical jet breakout time. Indeed, when we plot in fig. (|2J) the quantity dNgnn J dT^ 
instead of the traditionally shown dNoRs/d log(Tgo) (e.g. fig. HJ iKouveliotou et al.l Il993l ) 
this flat distribution is evident. The plateau appears over about an order of magnitude in 
duration around a few seconds, in the GRB duration distributions of BATSE, Swift and 
Fermi GBM, as depicted in fig. [2j The duration is characterized by T90 during which 90% 
of the fluence is accumulated. 

At the short end the distribution is rising towards shorter durations. This "bump" in the 
duration distribution is inconsistent with a Collapsar origin for most of the short duration 
GRBs. This simple conclusion is consistent with other evidence that a second, non-Collapsar, 



population of short duration GRBs exists with a different origin than the long 


er ones. (e.g. 


Kouveliotou et al. 


1993; 


Barthelmv et al. 


2005; 


Fox et al. 


2005 




Gehrels et al. 


2005; 


Nakar 



20071 ). 



To quantify the non-Collapsars' duration distribution we make joint fits to the overall 
duration distributions, including the Collapsar distribution at durations longer than the 
plateau. Although we are interested only in the short duration regime, where the duration 
distribution of Collapsars is flat, inclusion of the long end of the distribution is needed to 
determine the height of the plateau. To test the robustness of our result we fitted various 
functional forms for the distribution at long durations and verified that the height of the 
plateaus are consistent within the errors. The results presented here employ a plateau below 
the typical observed breakout time, Tg, and a power law with an exponential cutoff above 
it. For the non-Collapsars we find that the best fitted distribution function is a lognormal. 
Overall we fit the duration distributions to the function: 

1 T 90 <T B ^ 

3 -/3(T 90 -T s ) 



dN, 
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dTi 
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where the first term corresponds to non-Collapsars and the second one to Collapsars. 

We consider the data sets of BATSE 3, SwiftE, and Fermi GBM 0. We limit the data 
to the duration regime of 0-200 sec, which is enough to obtain a good constraint of the 
plateau hight. We verified that changing this range to 0-1000 sec has no significant effect 
on our results. We fit each sample with a distribution function according to eq. (pQ). After 
using the normalization that the integral of dNcRB/dTgQ over the duration range equals the 
number of observed GRBs, we are left with seven free parameters. We obtain good fits with 
X 2 per degrees of freedom (DOF) of 0.9, 1.3, 1.1 for the BATSE, Swift and Fermi GBM 
respectively. The corresponding parameters are given in table [1] and Fig. [2] depicts the 
resulting distribution functions and the data. We find plateaus that extend up to Tb ~ 20 
sec in the BATSE and Fermi GBM durations distributions, and up to Tr ~ 10 sec in Swift . 



This is consistent with our expectations from the Collapsar model (IBromberg et al.ll2011al ). 



Table 1: Best fit parameters 

Detector A^c A* ° Tb (s) a j3 

BATSE 545 -0.5 ±0.1 1.32 ± 0.07 25.5±^ 19.4^ -0.33 ± 0.2 0.019 ± 0.003 

Swift 42 -1.5 ±0.5 1.5 ±0.4 10.0 ± 2.3 7.9 ± 3.5 -0.3 ±0.2 0.01 ± 0.003 

Fermi GBM 128 -1.5 ±0.6 1.9 ±0.4 8.8 ± 0.2 18.2+^ -1.2 ±0.36 0.008 ± 0.001 



3. The non-Collapsars probability function 

The probability that a GRB with a given T90 is a non-Collapsar is given by the fraction 
of non-Collapsars within the observed GRBs at a given duration: 

1 (lnT 90 - M ) 2 / dN GRB 



f(T 9Q )=A NC — e-^^~ , (2) 

Igoay 2% V ai 90 



where dNcRB / 'dTg is given by eq. (DO). To estimate the errors in f^c we simulate, for each 
one of the samples, distributions of T go drawn randomly from the best fitted distribution 
function dNcRB/dT 90 . We then bin the simulated data sets, and repeat the process of 
parameter fitting using eq.((T]) to obtain f^c ■ We repeat this processes 1000 times and look 



2 http://swift. gsfc.nasa.gov/docs/swift/archivc/grb_tablc, from April 21, 1991 until August 17, 2000. 



3 |http://gammaray.msfc. nasa.go v/batsc/grb/catalog/c urrent /\ from December 17, 2004 until February 
20, 2012. 

4 http://heasarc. gsfc.nasa.gov/W3Browse/fermi/fcrmigbrst. html, from July 17, 2008 until July 9, 2010. 
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Fig. 2.— The T 90 distributions, dN/dT 90 , of BATSE (red), Swift (blue) and Fermi GBM 
(green) GRBs, binned into equally spaced logarithmic bins. Bins with less than 5 events 
are merged with their neighbors to get more reliable statistical errors. Note that the quan- 
tity dN/dT is depicted and not dN/d\og(T) as traditionally shown in such plots (e.g., 
Kouveliotou et al.lll993l ). The combined best fitted distribution functions of both Collap- 
sars and non-Collapsars are shown with solid lines. The dotted and dashed lines depict the 
distributions of non-Collapsars and Collapsars respectively in each data set. 
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for the ranges of f^c that encompasses 68% of the cases. Fig. [3] depicts /nc (^9o) for the 
BATSE, Swift and Fermi GBM samples. The solid lines depict f^c , calculated from the 
observed data, and the blue region describes the la error estimate. Table [T] lists the T 90 
values that correspond to some selected probabilities for the three detectors. 

These results clearly show that the choice of Tgo = 2 sec as a threshold to identify 
non-Collapsars is suitable for BATSE, and possibly also for Fermi GBM. A BATSE (Fermi 
GBM) burst with T 90 = 2 sec has a probability > 70% (> 40%) to be a non-Collapsar. 
However, the probability of a similar Swift burst to be a non-Collapsar is only 0.16 ± 0.14. 
It is most likely a Collapsar! The level of false identification, for a given duration threshold, 
T t h can be seen in Fig. H] that depicts the integrated fraction of Collapsars (out of the total 
number of GRBs) with duration T 90 ^ T^. The total number of Collapsars with duration 
T 90 < 2 sec in the Swift sample is estimated to be 19 ± 5 out of 53 GRBs. Thus, an arbitrary 
Swift sample selected with the T t h = 2 sec criterion contains about than 40% Collapsars that 
have been misclassified as non-Collapsars. This should be compared with about 10% and 
15% of misclassified Collapsars in the corresponding BATSE and Fermi samples (see fig. H]) 
with the same T t h- The criterion T f ^ = 2 sec for selecting non-Collapsars in Swift is simply 
very bad for most studies of non-Collapsars. 

Any single criteria that should distinguish according to the durations between Collapsars 
and non-Collapsars should be detector dependent. A longer T t h increases the size of the 
SGRBs sample (that are supposedly non-Collapsars) but it increases at the same time the 
number of misidentified Collapsars in the sample. A shorter threshold yields smaller but 
cleaner samples. The specific choice of T t h should be considered for each study, balancing 
the need of a large sample with the importance of purity. A reasonable choice that should be 
adequate to many studies is choosing the threshold probability as f^c — 0.5. This reconciles 
between these conflicting requirements and allows us to classify both Collapsars and non- 
Collapsars with a single criterion. Adopting this probability we find that the corresponding 
Tg threshold values are: T th = 0.8 ± 0.3 sec for Swift , T th = 1.7±° i sec for Fermi GBM 
and Tth = 3.1 ± 0.5 sec for BATSE. The total number of misclassified Collapsars in samples 
selected according to these criteria constitute about 20% of the Swift samples and ~ 14% of 
the BATSE and Fermi GBM samples (fig. 0}. 



3.1. The non-Collapsar probability as a function of duration and hardness 



As already mentioned short GRBs are harder on average than long ones (jKouveliotou et al 

19931 ). It is natural to expect that the ratio of non-Collapsars to Collapsars and the proba- 
bility function, f^c , should increase with the GRB hardness. Therefore the combination of 
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Fig. 3. — The fraction f^c of non-Collapsars as a function of T go for BATSE, Swift & Fermi 
GBM (from top to bottom). This fraction represents the probability that a GRB with an 
observed duration T 90 is a non-Collapsars. The shaded regions represent the 68% confidence 
range. The red vertical lines mark the values of T90 where f^c — 0.5 (See table |2]for numeric 
values of T 90 that correspond to some selected f^c values). 
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Fig. 4. — The integrated fraction of Collapsars with duration T 90 ^ T th , for BATSE, Swift & 
Fermi GBM (from top to bottom). This is the fraction of Collapsars that are misclassified as 
non-Collapsars with durations shorter that the threshold, Tth- The shaded regions represent 
the 68% confidence limits. 
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Table 2: The T90 (sec) that corresponds to some selected Jnc values in the three satellites 
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duration and hardness provides a stronger way to distinguish between non-Collapsars and 
Collapsars. To examine the duration-hardness probability we divide the samples into three 
hardness subgroups: soft, intermediate and hard, and preform the same analysis on each sub- 
groups: We fit a duration distribution function and calculate the probability function, f^c 
, and its la variance. We consider two hardness thresholds: The soft and intermediate sub- 
groups are separated by the average hardness of Collapsars which we estimate using GRBs 
with T90 > 20 sec, where the contribution of non-Collapsars in all samples is negligible. The 
intermediate and hard subgroups are separated by the average hardness of non-Collapsars 
which is estimated using GRBs with T 90 < 0.5 sec. The spectral hardness of different satellite 
samples is quantified differently for each detector, depending on the available information 
for each database. For BATSE we use the hardness ratio parameter, HR 32 , defined as the 
ratio between the photon counts in energy channel 3 (100 - 300 keV) and energy channel 2 
(50 - 100 keV). In Swift and Fermi GBM samples we use the powerlaw index (PL) of the 
observed spectrum obtained by fitting a single powerlaw in the energy range 15 — 150 keV 
{Swift ) or 10 - 2000 keV (Fermi). Note that in the Swift sample, only ~ 87% of the GRBs 
have a spectral fit to a single power-law. The spectrum of the other 13% is fitted with a 
powerlaw+exponential cutoff, and we omit these bursts from this analysis. 

If the distribution functions of Collapsars and non-Collapsars don't depend strongly on 
the spectral hardness, then varying the hardness threshold would only change the relative 
ratio of Collapsars to non-Collapsars, while the overall shape of the duration distribution 
functions of the two populations would remain unchanged. To examine this we fit each 
hardness subgroup with the same distribution function as in the full sample of the corre- 
sponding detector. We only rescale A^c and Ac according to the number of non-Collapsars 
and Collapsars in the subgroup relative to their number in the full sample. We evaluate the 
ratio of non-Collapsars by dividing the number of GRBs with T90 < 0.5 sec in the subgroup 
with their number in the full sample. The ratio of Collapsars is evaluated in a similar way 
using GRBs with T 90 > 20 sec. We find good fits in all hardness subgroups with x 2 /dof 
of ~ 0.8 — 1.8. This good fit indicates a weak dependency of the duration distributions of 
Collapsars and non-Collapsars on the hardness. 
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Figure [6] depicts the observed duration distributions of the three hardness subgroups 
of each detector. The harder subgroups have more prominent 'bumps' at short durations 
together with relatively lower plateaus that become visible only at longer durations. This 
is the expected behavior if the fraction of non-Collapsars increases with the GRB hardness. 
The dN/dTgo distributions and their ^/DOF values are shown in fig. [61 The resulting 
probability functions of each subgroup are shown in figs. [7]|9j In table [3] we list the T 90 values 
that correspond to specific values of f^c in eac h subgroup. In the hard BATSE subgroup 
non-Collapsars dominate the distributions up to T 90 = 7.81^ sec, while in the hard Swift 
and Fermi GBM samples non-Collapsars dominate up to Tgo = 2.8^}'o and T 9 o = 5. 4^0 
sec respectively. The transition between Collapsars and non-Collapsars in the intermediate 
subgroups roughly follow the same values as in the complete samples: Tgo = 2.9 ± 0.4, 
T 90 = 0.6±o;3 and T 90 = 1.6^:1 sec for BATSE, Swift and Fermi GBM respectively. In the 
soft subgroups non-Collapsars dominate only up to Tg = 1-IIqJ sec in BATSE, Tg = 0.31q; 2 
sec in Swift , and up to Tg = 0.6^o;3 sec in Fermi GBM. 
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Table 3: The Tc 
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at different f^c in the three hardness subgroups for each satellite. 
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The probability functions we obtained here can be used to classify the GRBs detected 
by BATSE, Swift and Fermi GBM according to their duration and hardness. For GRBs that 
cannot be assigned to one of those hardness subgroups (e.g the 13% of Swift GRBs whose 
spectra are fitted with a powerlaw+exponential cutoff) the classification can be done using 
the overall probability function of the complete Swift sample. For GRBs that have been 
detected by HETE and Integral we recall that the spectral window observed by HETE is 
8 — 500 keV, which is closer to the Swift /BAT while Integral, on the other hand, observes 
at a spectral range of 15 keV - 10 MeV, which is closer to BATSE's range. As a first 
approximation one can use the corresponding f^c values of these detectors. 

In Appendix |A] we collect all Swift GRBs with T 90 < 2 observed to date. The table in- 
cludes also GRBs with a hard short spike plus a soft extended emission and a number of other 
GRBs with duration > 2 sec that are sometimes considered as possible non-Collapsars. For 
each GRB we calculate the probability to be a non-Collapsar from the duration and power- 
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law index ,/jvc {Tgo, PL). For those GRBs with a spectral fit of a powerlaw+exponential 
cutoff, we calculate the probability to be a non-Collapsar from the duration alone, f^c (^90)- 
Important GRBs are emphasized with bold text. The table also includes a few important 
GRBs detected by HETE or Integral. For these GRBs we estimate f^c (T90) using the prob- 
ability functions of Swift and BATSE respectively. This table can be used to evaluate the 
contamination by Collapsars in present samples of SGRBs and to select low contamination 
samples in future studies. 



4. Consistency checks with studies of contaminated samples 



The commonly used criterion to distinguish Collapsars from non- Collapsars is the dura- 
tion (T90 ^ 2 sec). This criterion is applied to GRBs that are detected by all 7-ray satellites 
including Swift , that supplies the largest number of well localized short duration GRBs. As 
we have shown earlier Swift GRBs with T go > 0.8 sec have high probability to be Collap- 
sars. This could have led to a "Collapsar contamination" in current Swift samples of SGRBs 
that are based on the 2 sec criterion, and might have affected the r esults of studies based 
on these samples. Interestingly such studies (e.g. lBergerll2009l . l201ll ). have shown that the 
environments of SGRBs are different than the environments of LGRBs. Whereas LGRBs 
are associated with inten sive star formation, a ris e in low metallically irregular star form- 
ing galaxies (see however Levesque et al. 2010bl lcl; ISavaglio et al. 2012 , for examples of high 
metalicity LGRB hosts) and are concentrated towards star forming regions in their galaxies. 
SGRBs are associated with a broad distribution of galaxy types and arise in hosts with a 
broad range of star formation rate and metallicities and show a larger scatter in the distance 
distribution from their hosts' centers. One may wonder how these results are consistent with 
our claim that the 2 sec classification is not valid for the Swift sample. 



Table H] lists the GRBs and the host galaxy characteristics used in the iBergerl (120091 . 



201 ll ) sample. It also includes the probability that the associated SGRB are non- Collapsars 
(based the combination of duration and power-law index). In addition to eight 'classically 
selected' SGRBs (T 90 < 2 sec) this sample includes also four GRBs with T 90 > 2 sec. 
These GRBs are characterized by a short hard initial spike followed by a long tail of softer 
emission. They are often consid ered as non- Collapsars since their initial spikes resemble a 
classical SGRB (see lNakarll2007l ) . However, since the overall duration is not well defined our 
classification scheme cannot attribute a non-Collapsar probability to these bursts. 

Even though our probabilistic approach is incapable of determining whether a specific 
burst is a Collapsar or not, a clear picture emerges from table HI Four out of eight classifi- 
able bursts are non-Collapsars at very high probabilities. Two bursts are almost certainly 
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Collapsars while the last two are marginal: the probability of each one of those two to be a 
Collapsar is larger than 60%, however the probability that both are Collapsars is less than 
50%. These fractions are consistent with what is expected, according to our analysis, from 
a Swift sample with a 2 sec criteria for which ~ 60% of the bursts should be non- Collapsars 
and the rest Collapsars. 

Within the sub-sample of four non-Collapsars we observe a large spread in SFRs, in spe- 
cific SFRs and in galactic luminosities. Distances from the center of the host have typically 
large observational error, but at least one is quite far from the center (~ 44^3 kpc). There 
is not enough data to determine the metallicity. These results show a large spread in the 
observed quantities, in a large contrast with the rather narr owly distributed hos t properties 
of LGRBs (Collapsars). This is similar to the conclusion of iBergerl (120091 . 1201 ll ). It demon- 
strates that also when a less contaminated, but smaller, sample is examined non-Collapsars 
hosts have a different distribution than Collapsar hosts and consequently that the two pop- 
ulations have different progenitors. On the other hand, as expected, the properties of the 
hosts of the two Collapsar candidates are fully consistent with those of typical LGRB hosts. 
Finally, the properties of the hosts of the two bursts with marginal classification are also 
consistent with being either Collapsar or non-Collapsar hosts. 

The conclusion that the properties of the non-Collapsars' hosts are widely distributed 
whereas those of the Collapsars ' hosts a re na r rowly distributed implies that our classification 
is consistent with the results of lBergerl (120091 . l201ll ) even though the latter are based of a sig- 
nificantly contaminated sample. A wide distribution contaminated by a narrowly distributed 
population retains it basic feature of a wide di stribution, and this is what happens here. The 
non-Collapsars within the IBergerl (120091 . 1201 ll ) SGRB sample are numerous enough to result 
in a wide distribution that is significantly different from the one of Co l lapsar s. However, 
while our conclusions are in line with the basic results of IBergerl (120091 . 1201 ll ). the details 
of the distribution, such as the ratio of high SFR to low SFR hosts or the distribution of 
distances from center, are influenced by the contamination and a quantitative study of the 
distribution of the host properties should take this factor into account. 

The possible effects of contaminating Collapsars on studies of properties of SGRBs vary 
from one study to another. Different samples have different contaminations and different 
properties are influenced differently. The probabilities given in appendix A can be used 
to evaluate the likelihood that different bursts are non-Collapsars or Collapsars and with 
these to estimate the quality of a specific sample and the significance of results based on 
this sample. In general one should proceed with care before adopting simply the results of 
a study of an SGRB sample as reflecting the properties of non-Collapsars. In this context 
it is interesting to mention a few GRBs that play a major role in the current view of non- 
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Collapsar properties. GRB 060121 and GRB 090426 are two SGRBs with a secure host at 
redshift ^ 2 that have led to the suggestion of a high redshift non-Collapsar population. GRB 
100424A has a redshift of z — 1.288. All other SGRBs with secure redshift are at z < 1. We 
find that the probabilities that these bursts are non-Collapsars are 0.171^15; 0.101°;^ and 
0.08ig;^ for GRBs 060121, 090426 and 100424A respectively. Surely, one cannot establish 
a new population of high redshift non-Collapsars based on these events. Another pivotal 
burst is 051221A; the only SGRB to date with a clea r simultaneous optical/X-ray break in 



its af terglow, which is used to measure its beaming (ISoderberg et al. 



2006; Burrows et al. 



+0.08 Thig 



20061 ). We find that the probability that GRB 051221A is a non-Collapsar is 0.18_ ( 
highlights our ignorance of the collimation (if there is any) of non-Collapsar outflows. It also 
highlights the fact that no firm conclusion can be drawn on non-Collapsars based on a single 
burst that is classified using its high energy emission properties alone. 



Table 4: The sample of [Berger] (|2009|, l201l|) SGRBs 
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a Swift GRBs with a single power-law spectral fit are assigned a probability f^c (T 90 ,PL). 

Other GRBs can only be assigned a probability Jnc (Tbo). 
b A GRB detected by HETE, Jmc (T90) is estimated using Swift probability function, 
f GRB with an extended softer emission 

* The metallicity is measured by the ratio of Oxygen to Hydrogen lines. The range of values 

of 8.2 — 8.9 shown in the table corresponds to ~ 0.3 — 1.6-Z©. 
R eferences: 1) iBergerl |2009[); 2) libxet al.lj2005j) ; 3) lFong. Berger. fc Foxl |2010l) ; 
4) iBerger et all |2009h : 5) ISoderberg etall |2006l ) 
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5. Summary 

GRBs are widely classified as long and short, according to their duration T go ^ 2 sec, 
based on the general belief that this observational classification is associated with a physical 
one and that the two populations have different origins: long GRBs are Collapsars and short 
ones are non- Collapsars (possibly arising from neutron star mergers, but at present, this 
association is still uncertain). This classification scheme is known to be imperfect due to the 
large overlap in the duration distribution between the two populations. It is also used for 
all de tectors, alth ough it is known that any classification scheme depends on the detector 



(e.g., iNakarl 120071 ). The problem with this method is that, first it is impossible to know 
how trustable are results that are based on a single classified event. Second, the level of 
contamination in any studied sample is unknown. The main reason for this flawed practice 
is simply the lack of a reliable and quantifiable classification scheme. This is what we provide 
i n this paper. Based o n a physically motivated model we have shown in an earlier study 



(IBromberg et al.ll2011al ) that at short durations the Collapsar distribution is flat, up to a 
typical duration of ~ 20 sec. This enables us to recover the non-Collapsar distribution from 
the overall duration distribution and to assign probability that a burst with a given duration 
and hardness is a non-Collapsar. 

We carry out this analysis for three major GRB satellites, BATSE, Fermi and Swift . 
We first find the probability that a burst is a non-Collapsar based on its duration alone, f^c 
(Tgo). We find that it depends strongly on the observing satellite and in particular on its 
spectral window. For a given duration the probability that a BATSE burst is a non-Collapsar 
is larger than the probability that a Swift burst is a non-Collapsar. A useful threshold 
duration that separates Collapsars from non-Collapsars is that where f^c (^9o) = 0.5. We 
find that it is T 90 = 3.1±0.5 sec in BATSE, T 90 = 1.7i{j;| sec in Fermi GBM and T 90 = 0.8±0.3 
sec in Swift. 



As short GRBs are harder on average than long ones (jKouveliotou et al.l 119931 ). it is 
natural to expect that GRBs with a hard spectrum have a higher probability to be non- 
Collapsars than softer ones. Thus, a better classification can be achieved by considering 
the hardness, in addition to the duration. We separate the sample of each satellite to three 
sub-samples based on the bursts hardness and repeat the analysis. Not surprisingly there 
are fewer non-Collapsars in the soft subgroups and more in the harder ones. Interestingly 
the duration distributions of both Collapsars and non-Collapsars depend only weakly on the 
hardness and only the relative normalization between the two groups varies as we consider 
subgroups of different hardness. As there are more non-Collapsars in the hard subgroups, 
non-Collapsars dominate in these subgroups even at relatively long durations. For example 
In the hard BATSE subgroup the probability, f^c , that a burst is a non-Collapsar remains 
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> 0.5 up to durations T go ~ 8 sec. In the hard Fermi GBM and Swift subgroups f^c > 0.5 
up to ~ 5 and ~ 3 sec respectively. A soft GRB, on the other hand, is more likely to 
be a Collapsar. In this case f^c > 0.5 up to T 90 ~ 1 sec for BATSE's soft subgroup, up 
to Tg — 0.6 sec in Fermi GBM and only up to T 90 ~ 0.3 in Swift 's subgroups. These 
values should replace the average values as dividing durations between Collapsars and non- 
Collapsars, whenever hardness information is available. In particular for Swift , 2.8l^o, 
0-6±8:i, and 0.3 

-o.2 sec should replace the value of 0.8 ± 0.3 sec for the hard, intermediate 
and soft subgroups respectively. Our results well agree with the overall behavior seen when 
comparing different satellites. Swift 's window is much softer than BATSE's and Fermi's 
and the transition in Swift 's overall sample between non-Collapsars and Collapsars occurs 
at shorter durations relative to the other satellites. This is a general pattern seen in both 
the overall sample and in the hardness subgroups. 

We find that the transition between Collapsars and non-Collapsars is not sharp and that 
there is a large overlapping region where both Collapsars and non-Collapsars co-exist. There 
are short durations Collapsars with durations shorter than 1 sec as well as non-Collapsars 
at observed durations as long as 10 sec. The traditional method to divide bursts to "long" 
and "short" according to a sharp observed duration criteria: T 90 ^ 2, introduces both 
"false positive" and "false negatives" when we interpret duration as a proxy for a different 
physical origin. The choice of the division criteria should depend on the detector's observing 
windows but it should also depend on our tolerance for contamination by "falsely" classified 
bursts. When interested in Collapsars, the solution is trivial. Choosing a conservative 
large duration will eliminate a few short duration Collapsars but will results in a sample 
containing practically only Collapsars. The small number of short duration bursts makes it 
difficult to adopt a similar conservative policy for them and the classification criterion should 
be chosen carefully in each study. Finally, our results show clearly that no high significance 
result concerning non-Collapsars can be derived based on a single burst, which is classified 
according to its high energy properties alone. 

Next we examine the implication of the currently used criterion Tgo ^ 2 on the different 
satellite samples. It is conservative for BATSE, where Only 10% of bursts that are shorter 
than 2 sec are Collapsars. One can consider a BATSE (Tg < 2 sec) sample as reasonably free 
of false positives. The corresponding fraction of "false positives" for Fermi is higher (20%) 
but still acceptable for many purposes. However this criteria is not good for Swift . About 
40% of Swift bursts with Tgo < 2 sec, that have been traditionally classified and studied as 
SGRBs are Collapsars. Thus the standard and commonly use sample of Swift GRBs with 
Tg ^ 2 sec which is the source of the only sample of well localized short GRBs is heavily 
contaminated with Collapsars! This must have influenced the results of non-Collapsars 
studies that are based on Swift GRBs. Interestingly, this Collapsar contamination didn't 
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affect qualitatively the main conclusion concerning non-Collapsar hosts (lBergerll2009l . 120111 ) . 
namely the observation that these hosts have a wide distributions of SFR, luminosities and 
metallicities and that the conclusion that the positions of non-Collapsars has a wide spread 
within the host galaxy. Such distributions are significantly different than those of Collapsar's 
host. However, quantitative features of these distributions must have been distorted. 

While the complete implications of our results on studies of non-Collapsars is beyond 
the scope of this work, there are three important points that stand out. (i) There is no 
convincing evidence for high redshift non-Collapsars. All the bursts with secure redshift that 
are non-Collapsars at high probability are at z < 1. (ii) There is no convincing evidence for 
beaming in non-Collapsars. GRB 051221A is the only SGRB that show a mult i- wavelength 
afterglow break, that is interpreted as a jet break and is considered as the strongest evidence 
for beaming in non-Collapsars. However, our results show that the probability that this 
burst is indeed a non-Collapsar is only 0.18^oii. Apparently, non-Collapsars may or may 
not be beamed as far as we currently know, (iii) The duration of a non-negligible fraction 
of the non-Collapsars is 10 s and even longer. This implies (under most GRB models) that 
the central engine of these events works continuously in the mode that produces the initial 
hard GRB emission for that long. This fact should be accommodated by any model of 
non-Collapsar central engine. 

This research is supported by an ERC advanced research grant, by the Israeli Center for 
Excellence for High Energy AstroPhysics (T.P.), by ERC and IRG grants, and a Packard, 
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a Swift GRBs with a single power-law spectral fit are assigned a probability Jnc (T90, PL) 



Other GRBs can only be assigned a probability f^c (Tgo)- 
b A GRB detected by HETE, f^c (T90) is estimated using the Swift probability function 
c A GRB detected by Integral, fjyc (T90) is estimates using the BATSE probability function 
d The spectral fit of the 7-ray photons is a power-law with an exponential cutoff, 

Inc (T$q,PL) cannot be calculated for this burst and f^c (T90) is used instead. 
' A GRB with an extended soft emission, no fxc is assigned. 

* Unsccure redshift, based on an association of a galaxy within the XRT error circle. 

Redshift references: fl jProchaska et ail j2005ah ; [Gehrels et al.l ||2005|) : (2 lVillasenor et alj j2005h ; lFox et ail j2005h ; 
O lBerger et ail [120051 1 ; IProchaska et all [|2005bh ; f4 lGehrels et ail J2005h ; iBergeil J2005h : [Foley, Bloom, fc Chenl l l2005h ; 
(5 lBerger fc Soderberd J2005h ; ISoderberg et all lb006l l; (6~lde~ Ugarte Postigo et al, ] d2006h ; I Levan et al. I ll2006l) ; 
(7 lBloom et ail j2006h : (8 lOfek et al.l boodl ; iLevesque fc Kewlevl J2007h ; O lPrice. Berger. fc FoJ feOOCft ; iFugazza et al.l j2006h : 
(lO lCucchiara. Cannizzo. fc BergeJ j2006h : (ll lBergeJ <2007ah ; (12 lBergeJ <2006al . l2007bh ; (13 lCenko et ail d2006h : 
(14 lBergeJ fe006bh ; (15 lPerlev et al.l J2007h ; 16 lGraham et al.l j2007h ; (17 lCucchiara et al.l j2007h ; ICovino et ail J2007f) ; 
(18 lD'Avanzo et ail i2007l) ; [Berger. Morrell. fc Rothl <2007ch ; (19 lRowlinson et ail feoid) ; (20 lLevesQue et ail (l2009h ; 
iThoene et ail j200Sh ; (2l lRau. McBreen. fc KruehleJ J2009h ; (22 lFong et al.l boill) ; (23 lThoene et ail bold) 




Fig. 5. — The hardness ratio of BATSE GRBs (top left) and the powerlaw index of Swift 
GRBs (top right) and Fermi GRBs (bottom) as a function of Tgo- 



-26 - 




Fig. 6. — dN/dT 90 , of the three hardness subgroups of the BATSE (upper left) Swift (upper 
right) and Fermi GBM (lower) samples, binned into equally spaced logarithmic bins. Bins 
with less than 5 events are merged with their neighbors to reduce statistical errors. The best 
fitted joined distribution functions are marked with solid lines. 
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Fig. 7. — The fraction, f^c , of non-Collapsars as a function of the observed duration, T 90 , 
in the 3 hardness subgroups of BATSE. 
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Fig. 8. — Same as fig. [7] for Swift 
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Fig. 9. — Same as fig. [7] for Fermi GBM 



